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ABSTRACT 

High-resolution, multi-wavelength, and time-domain observations of the Galactic centre black 
hole candidate, Sgr A*, allow for a direct test of contemporary accretion theory. To date, all 
models have assumed alignment between the accretion disc and black hole angular momen- 
tum axes, but this is unjustified for geometrically thick accretion flows like that onto Sgr A*. 
Instead, we calculate images and spectra from a set of numerical simulations of accretion 
flows misaligned ("tilted") by 15° from the black hole spin axis and compare them with mil- 
limetre (mm) to near-infrared (NIR) observations. Non-axisymmetric standing shocks from 
eccentric fluid orbits dominate the emission, leading to a wide range of possible image mor- 
phologies. These effects invalidate previous parameter estimates from model fitting, including 
estimates of the dimensionless black hole spin, except possibly for very low values of spin or 
tilt (upper limits of a < 0.3 or /? < 15°). At 1.3mm, the black hole images still have crescent 
morphologies, and the black hole shadow may still be accessible to future very long base- 
line interferometry (mm-VLBI) observations. Shock heating leads to multiple populations of 
electrons, some at high energies (Te > lO^^K). These electrons can naturally produce the 
observed NIR flux, spectral index, and rapid variability ("flaring"). This NIR emission is un- 
correlated with that in the mm, which also agrees with observations. 

These are the first models to self-consistently explain the time-variable mm to NIR emis- 
sion of Sgr A*. Predictions of the model include significant structural changes observable with 
mm-VLBI on both the dynamical (hour) and Lense-Thirring precession (day-year) timescales; 
and ~ 30 — 50/ias changes in centroid position from extreme gravitational lensing events dur- 
ing NIR flares, detectable with the future VLT instrument GRAVITY. We further predict that 
multi-wavelength monitoring should find no significant correlations between mm and NIR/X- 
ray light curves. The weak correlations reported to date are shown to be consistent with our 
model, where they are artifacts of the short light curve durations. If the observed NIR emis- 
sion is caused by shock heating in a tilted accretion disc, then the Galactic centre black hole 
has a positive, non-zero spin parameter (a > 0). 

Key words: accretion, accretion discs — black hole physics — radiative transfer — relativity 
— galaxy: centre 



1 INTRODUCTION 

Due to its proximity, Sgr A* is the most intensively studied su- 
permassive black hole candidate and the largest in angular size. 
Recent ver y long baseline interferometry observat ions at 1.3mm 
(mm-VLBI [Poeleman et"ani2008l : lFish et alj|201lh have detected 
source structure on event-horizon scales (tens of microarc seconds, 
Has). Future measurements with this Event Horizon Telescope may 
detect the black hole shadow, providing the first direct evidence 
for an event horizon ( lBardeen|[l973l Ipalcke. Melia & Agolll2000l : 
lDexteretal.ll2010h . 



The spectral energy distribution of Sgr A* rises from the radio 
to a peak in the submillimetre (mm "bump"). The radio emission 
is synchrotron radiation from non-ther mal electrons far from the 
black hole, either in the accretion flow (Yuan. Ouataert & N aravanI 
|2003|) or in a mildly relativistic jet (Falcke & Markoff 2000). The 
mm bump is well described by synchrotron radiation from hot (~ 
■|^gio-ii j^-j^ mostly thermal elec trons in the immediate vicinity of 
the black hole (r 5 - lOlvTHYjan Ouataert & NaravanI 120031: 
iMoscibrodzka et al.ll2009l:lDexter et al.ll2010h . 

In addition to 30 — 50 per cent variability in the mm 
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I Zhao et alj l2003l: lYusef-Zadeh et alj I2OO6I : lEckart et alj I2OO6I : 



iMarrone et al. 2008^, Sgr A* i s observed to "fl are" rapidly 
in t he NIRloenzel et ah l2003l : lOhez et alj l2004h and the X- 
ray ( iBaganoff et alj l200lh with order of mag nitude amplitudes . 
The NIR/X-ray flares occur simultaneously dEckart et al.l l2004l) . 
and weak correlations in simultaneous light curves have sug- 
ested possible lags be tween NIR/X-ray and mm variability 
Yusef-Zadeh et alj|2008h . The NIR (X-ray) emission is either syn- 
chrotron radiation from energetic ele ctrons, or inverse C ompton 
scattering of mm (NIR) seed photons jMarkoff etaT]|200lh . 

No theoretical model of Sgr A* self-consistently accounts for 
the observed time-variable multi- wavelength emission. Radiatively 
ineffeci ent accretion flow (RIAFmian. Ouataert & Naravan 2003.) 
and jet dFalcke & MarkofjbOOd) models can describe the multi- 
wavelength "quiescent" spectrum. RIAF models also provide ex- 
cellent fits to existing mm-VLBI data iBroderick et al. 2009, 201 ll). 

Models for the NIR flares include coherently orbiting in- 
homogeneities ("hotspots" IBroderick & Loebl l2006h in the inner 
radii of the accretion flow. Recently, correlated multiwavelength 
flares have instead favored an adiabatically expanding blob model 
(Yusef-Zadeh et al. 2009; Eckart et al. 2012), which can explain the 
reported time lags between the NIR/X-ray and mm flares. The blobs 
could be formed from non-th ermal events such as ma gnetic recon- 
nection in the accretion flow jPodds-Eden et al]|20ICll) . 

All of these models are non-relativistic. They neglect 
the magnetic fields responsible for angular momentum trans- 
port and accretion via the magnetorotational instability (MRI, 
iBalbus & HawievI Il99lh . Magnetohydrodynamic (MHD) simu- 
lations can provide a more physical description of the accretion 
flow. They have been used to model the synchrotron emis- 
sion from Sgr A*, either in three spatial dimensions with a 



pseudo-Newtonian potential jOhsuga. Kato & Mineshigd 20051 : 
Goldston. Ouataert & IgumenshchevI 20051 ; Huang etaU 20091 ; 



Chanetal. 2009), or in fu ll general relativity (GRM HD) 
in two dNoble etal] |2007|; iMoscibrodzka et al.l |2009|) or 
three jPexter. Agol & Fragilel l2009t JPexter et al.l [2OIOI ; 
IShcherbakov. Penna & McKinnevI l20ld ; IPolence et all I2OI2I) 
dimensions. Non-relativistic simulations are especially inappropri- 
ate for modeling the mm/NIR/X-ray emission, which originates 
in the innermost portion of the accretion flow where relativistic 
effects are strongest. Axisymmetric simulations cannot sustain the 
MRI, and cannot accurately model variability. 

Existing models from 3D simulations provide an ex- 
cellent de scription of mm observations including polarization 
(IShcherba kov. Penna & McKinnevI 201(]|), variability an d mm- 
VLBI (iDexter. Agol & Fragilel l2009l ; IPexter et al.l l2010h. These 
mode ls also produce large amplitude NIR flares i Dolence et alj 
|20I2|) . but with much lower fluxes and steeper spectral indices than 
observed. Although current global GRMH D simulations are factors 
of a f ew in resolution from convergence dHawlev. Guan & KrolikI 
l20Ill) . the r esulting spectra of Sgr A* are fairly insensitive to reso- 
lution l lShi okawa et al. 2012). 

There are two major limitations to current models. First, the 
electron distribution function is highly uncertain. Most theoretical 
models adopt a thermal distribution for the mm bump, and some 
add a non-thermal power law tail. Models from simulations use a 
thermal distribution with a constant electron-ion temperature ratio, 
Ti/Te, allowing the electron temperature to be calculated from the 
output density and pressur e ( Goldston. Ouataert & IgumenshchevI 



I2OO5I ; IMoscibrodzka et al.ll2009l) . Simple prescriptions for allow- 
ing this ratio to vary spatially lead to only minor changes in the 
images/spectra. This parameter can be calibrated either using phys- 



ically motivated sub-grid models for electron heating dSharma et alj 
I2OO7I) or through direct particle-in-cell calculations of MRI turbu- 
lence (|Riquelme et al. 2012). The second major limitation to all ex- 
isting theoretical models of Sgr A*, and the subject of this article, 
is the assumption that the accretion flow angular momentum axis is 
aligned with the black hole spin axis. 

Geometrically thick accretion flows like that onto Sgr A* do 
not align with a central bla ck hole as geometrically thin discs do 
dBardeen & Pettersonlll975l). S^r A* is c urrently fed by the winds 
of young, massive stars ( Ouataenll2004t ). whose net angular mo- 
mentum axis is not expected to be aligned with the spin axis of Sgr 
A*. The only way for alignment to occur is then for the black hole 
to accrete an amount of material with constant angular momentum 
orientation comparable to its mass. For this to occur in a Hubble 
time requires an accretion rate ~ 10~'*Myr"^, an order of mag- 
nitude larger than the estimated a ccretion rate at the Bondi radius 
dYuan. Ouataert & Naravanll2003l) and several orders of magnitude 
larger than that onto Sgr A* (10~^ - 10~^M(7 )yr~ HAitken et all 
I2OOOI ; iBower et alj |20(}3| ; iMarrone et J] 120071) . Therefore, it is 
highly likely that the accretion flow is misaligned ("tilted"). 

The expected disc misalignment has extraordinary conse- 
quences for the dynamics and observable properties of a black 
hole accretion flow. Lense-Thirring precessi on leads to tilts and 
warps that are oscillatory functions of radius dlvanov & lUarionovl 
and gravitational torques cause the en tire disc to pro- 
cess essentially as a solid body dFragile et al] [20073. Eccentric 
orbits in the disc converge, leadin g to non-axisymmetric stand- 
ing shocks dFragile & BlaesI 1200 Sl) . Excess angular momentum 
transport from the shocks truncates the accretion flow outside o f 
the marginally stable orbit ('Fragile"2009'; Dexter & Fr agilell201ll) . 
The asymmetry of the accretion flow leads to images and spec- 
tra that depend strongly on the observer azimuthal viewing angle 
dDexter & Fragilell20I II) . 

In this work we calculate time-dependent images and spectra 
from the only published GRMHD simulations of tilted black hole 
accretion discs dFragile et al. 20071 l2009tlFragile|[2009) , and com- 
pare them with mm to NIR observations of Sgr A*. We review the 
physics of tilted accretion discs in ^and describe the numerical 
simulations in Sj3] The methods used to construct radiative models 
from the simulations are outlined in ij4| These radiative models are 
then fit to existing mm-VLBI and spectral observations in [J5] and 
the mm images and mm to NIR spectral and variable properties 
of the best-fitting models are analyzed in !j6| The implications and 
limitations of this study are discussed in ^ and the major results 
are summarized in Sj8] 



2 TILTED BLACK HOLE ACCRETION DISCS 

Whenever a black hole accretion disc is tilted, it experiences dif- 
ferential Lense-Thirring precession owing to the frame dragging of 
the rotating black hole. In a disc, Lense-Thirring precession is able 
to build up over many orbital periods, so it can affect the structure 
well beyond radii normally associated with relativistic effects. 

The effect of Lense-Thirring precession in a disc is to cause it 
to warp, and this warping can propagate through the disc in either 
a diffusive or wave-like manner. In the diffusive case, the warping 
is limited by "viscous" responses within the disc, and the Lense- 
Thir ring precession dominates out to a unique transit ion radius 
rBP dBardeen & Pettersonll975l ; lKumar & Pringlell985l) . inside of 
which the disc is expected to be flat and aligned with the black-hole 
midplane, and outside of which it is also expected to be flat but in 
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Figure 1. Images as a function of Ti/Te for the 90h, 315h, and 515h simulations. The colours are scaled linearly with a dynamic range of 60, and the images 
are 160 X 160/ias in all cases. Images from untilted simulations such as 90h become optically thick and emit from a growing photosphere with increasing 
Ti/Te. Images from tilted simulations do not become uniformly optically thick, and can even become smaller at large Ti/Te, as in the case of 515h, due to 
the presence of multiple electron populations with varying electron temperatures at distinct spatial locations in the accretion flow. 



a plane determined by the angular momentum of the gas reservoir. 
This "Bardeen-Petterson" configuration is expected for Keplerian 
discs whenever the dimensionless stress parameter a is larger than 
the ratio H /r, where H is the disc semi-thickness. Given that a 
is usually considered to be significantly less than one, this requires 
geometrically thin discs. 

In thicker discs, such as would be expected in low-luminosity 
sources such as Sgr A*, warps propagate in a wave-like manner, 
and the disc does not go into the Bardeen-Petterson configuration. 
Instead, the tilt of the disc can be an oscillatory function of radius, 
with the amplit ude of the oscillations gro wing as one approaches 
the black hole jivanov & Illarionovlll997h . These oscillations are 
stationary, however, so they cannot act to absorb the torque of the 
black hole. For thick discs, this torque instead causes the disc to 
precess as if it were a solid body. 

The warping of tilted thick discs creates unbalanced pressure 
gradients that can drive significant latitude-dependent epicyclic 
motion within the disc. The induced motion of the gas can be co- 
herent over the entire scale of the disc, with magnitudes that are 
substantial fractions of the orbital velocity. Radial variations in the 
eccentricity of the associated fluid-element trajectories can lead to 
a crowding of orbit trajectories at certain locations within the disc. 
This results in local density enhancements akin to compressions. 
These compressions can be sufficiently strong to produce standing 



Table 1. Simulation Parameters 



Simulation 


a 


P 


Reference 


315h 


0.3 


15° 


Fragile (2009) 


515h 


0.5 


15° 


Fragile et al. (200^ 


715h 


0.7 


15° 


Fragile (2009) 


90h 


0.9 


0° 


Fragile et al. (2007) 


915h 


0.9 


15° 


Fragile et al. (2007) 



shocks within the disc, particularly in the vicinity of the black hole 
jpragile & BlaesllOOSt) . 

Tilted accretion disc models can be constructed semi- 
analytically usi ng height-integrated equati o ns in the so-called 
twiste d frame (l lvanov & Illarionovl Il997t IZhuravlev & Ivanovl 
I2OI ll) . However, these models lack the magnetic fields that lead 
to angular momentum transport in the disc and give rise to the ob- 
served synchrotron emission from Sgr A*. They also have an un- 
constrained vertical structure, leading to additional free parameters. 
Finally, they are stationary, and so cannot be used to study variabil- 
ity self-consistently. Instead, we employ 3D GRMHD simulations 
of misaligned black hole accretion discs. Simulations have fewer 
free parameters and higher physical fidelity. Few exist, however, 
preventing anything like a thorough exploration of the parameter 
space in dimensionless black hole spin, disc tilt, and geometric disc 
thickness. 
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Figure 2. Normalized probability distributions as functions of black hole spin (top left), observer inclination (top right), sky orientation (bottom left), and 
observer azimuth (bottom right) for each simulation separately and combined (solid lines). Unlike in aligned models ifiroderick et al.l2009ll201 lllDexter et alj 
lioiQ) the parameters of the viewing geometry are essentially unconstrained. 
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Figure 3. NoiTnalized probability density as a function of electron temperature (left) and accretion rate (right), marginalized over the other parameters. The 
electron temperature is taken to be the median in the region of largest emissivity for each model. The accretion rate is time-averaged for each model. In both 
cases, probability densities for models falling in each bin are added to make the histogram. The dotted histogram s in both panels are the result of assigning 
equal probability to each model. The viable range of accretion rates is much larger than in aligned simulations ( Dexte r. Agol & Fragilell200^ iDexter et al.l 
while the electron temperature remains well constrained to a narrow region around the observed brightness temperature. 



3 SIMULATIONS 

The numerical data used in our analys i s were taken from simu- 
lations presen ted in gile et al ] (l2007h . [Fragile et al.l (l2009l) . and 
[pragiie ( 2009). All of the simulations used the Cosmos-l-l- GRMHD 
code tAnninos. Fragile & Salmonson ,2005.) . Cosmos-l~l- includes 



several schemes for solving the GRMHD equations; in these sim- 
ulations, the internal-energy, artificial viscosity formulation was 
used. The magnetic fields were evolved in an advection-split form, 
while using a hyperbolic divergence cleanser to maintain an ap- 
proximately divergence-free magnetic field. The GRMHD equa- 
tions were evolved in a "tilted" Kerr-Schild polar coordinate system 
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{t, r, t}, (p). This coordinate system is related to the usual (untilted) 
Kerr-Schild coordinates {t, r, 6, (j)) thro ugh a sim ple rotation about 
the y-axis by an angle /3, as described in lFragile & Anninos (2005). 

The simulations were carried out on a spherical polar mesh 
with nested resolution layers. The base grid contained 32'^ mesh 
zones and covered the full An steradians. Two levels of refinement 
were added on top of the base layer; each refinement level doubling 
the resolution relative to the previous layer, thus achieving peak 
resolutions equivalent to a 128"^ simulation. In the radial direction 
a logarithmic coordinate of the form rj = 1.0 + ln(r/rBH) was 
used, where tbh is the black-hole horizon radius. In the angular 
direction, in addition to the nested grids, a concentrated latitude 
coordinate X2 of the form il) = X2 + ^{1 — h) sin(2x2) was used 
with /i = 0.5 to concentrate resolution toward the midplane of the 
disc. 

The simulations started from the analytic solution for an ax - 
isymmetric torus around a rotating black hole (IChakrabartilll985h . 
The inner radius of the torus was rin = 15M; the radius of the 
initial pressure maximum was rccntrc ~ 25M; and the power-law 
exponent used in defining the initial specific angular momentum 
distribution was q = 1.68. An adiabatic equation of state was as- 
sumed, with r = 5/3. The torus was seeded with a weak dipole 
magnetic field in the form of poloidal loops along the isobaric con- 
tours within the torus. The field was normalized such that initially 
/3mag = P/Pb ^ /3mag,o ~ 10 throughout the torus. The black 
hole was inclined by an angle /? = 0° or 15° relative to the disc 
(and the grid). From this starting point, the simulations were al- 
lowed to evolve for a time equivalent to 10 orbits at the initial pres- 
sure maximum, rcontrc, corresponding to hundreds of orbits at the 
ISCO. Table[T]summarizes the simulation parameters. 



4 RADIATIVE MODELING 

A radiative model is specified by three components: a dynamical 
model for fluid variables as functions of space and time; a model 
for the electron distribution function; and an emission model. In 
this work, the GRMHD simulations described above are used 
as the dynamical model. As in previ ous radiative models of Sgr 
A* from simulations (Goldston, Ou ataert & IgumenshchevI boOSi : 
iMoscibrodzka et alj2009t ; iDexter et alj2010l) . we emplov a thermal 
Maxwell- Juttner distribution for the electrons, assuming a constant 
ion-electron temperature ratio, Ti/Te, as a crude approximation to 
the electron distribution function in a collisionless plasma. This ra- 
tio is a free parameter in the models. Finally, we assume the emis- 
sion is entirely from sync hrotron radiation, and use the unp olarized 
emission coefficient from lLeung. Gammie & NobldfeOlU) . 

With these assumptions, observables can be calculated from 
the simulations using ray tracing. Starting from an observer's cam- 
era, rays are traced backwards in time toward the black hole (as- 
suming they are null geodesies) using the public code GEOKERR 
(Dexter & Agol 2009). In the region where rays intersect the ac- 
cretion flow, the radiative transfer equation is solved along the 
geodesic using the code GRTRANS tDexteill20Ilh . which then rep- 
resents a pixel of the image. This procedure is repeated for many 
rays to produce an image, and at many time steps of the simula- 
tion to produce time-dependent images (movies). Light curves are 
computed by integrating over the individual images. Repeating the 
procedure over observed wavelengths gives a time-dependent spec- 
trum. 

To calculate fluid properties at each point on a ray, the space- 
time coordinates of the geodesic are transformed from Boyer- 



Lindquist to the tilted Kerr-Schild coordinates used in the simula- 
tions. Since the accretion flow is dynamic, light travel time delays 
along the geodesic are taken into account. Data from the sixteen 
nearest zone centres (eight on the simulation grid over two time 
steps) were interpolated to each point on the geodesic. The sim- 
ulations are scale free, and must be scaled to cgs units in order 
to calculate emission and absorption coefficients. Fixing the black 
hole mass (Afbh ~ 4 x 10*^ Mq for Sgr A*) sets the length and 
timescales. The total mass in the accretion flow is a free parame- 
ter, and its choice is equivalent to fixing the time-averaged accre- 
tion rate. To convert image intensities to flux, we adopt a distance 
D — 8kpc to the Galactic centre. 

For each simulation, images are produced over a grid of pa- 
rameters at 0.4mm and 1.3mm: dimensionless black hole spin, a 
(or equivalently simulation); time-averaged accretion rate, AI; ion- 
electron temperature ratio, Ti jT^; observer inclination, i, azimuth, 
and sky orientation (position angle measured E of N), ^. These 
images (including the effects of interstellar scattering, iBower et al 
2006) are fit jointly to mm-VLBI observations (Do eleman et al 
i2008i ; lFish etal.l201lh and spectral measurements (Ma rronel2006t) . 
This gives the probability of observing the measured values from a 
particular model, which is converted to the probability distribution 
as a function of the mo del paramet ers using the method described in 
lBrodericket"aLl ( l2009h and used in lDexteretal.l ( l2OI0l) . This proce- 
dure gives a joint probability distribution as a function of the above 
parameters. To estimate individual parameters, we marginalize over 
the others using uniform priors on a, Ti/Te, (/>o, and ^; a logarith- 
mic prior on Af , and a sin i prior on the inclination angle. 

In practice, the range of M is used to bracket the range of ob- 
served mm (230 GHz) fluxes from Sgr A*: ~ 2 — 5Jy. Achieving 
convergence requires a fine enough grid so that multiple images 
have mm fluxes close to those observed by mm-VLBI (~ 2.4Jy) 
and spectral observations (~ 3.5— 4.25Jy). In untilted simulations, 
this is possible with a fairly coarse grid of 10 M values, since the 
image structure is nearly unchanged as the flux varies by 30 — 50 
per cent (see [j6jT). Variability then acts to sample the total flux 
more finely. This is not the case in tilted simulations, where the 
image structure changes drastically on short timescales. In proba- 
bility distributions over both observer time and position angle (^), 
reasonable convergence is obtained when using 20 values of M. 

The observer azimuth, 00, is a new parameter for tilted sim- 
ulations, since the accretion flow is no longer azimuthally sym- 
metric on times much longer than the orbital time. Instead, non- 
axisymmetric structure in these simulations averages out only on 
the much longer precession timescale0 Also unlike untilted sim- 
ulations, th e non-axisymmetric stru cture dominates images from 
tilted discs dDexter & Fragilell20Ill) . We find that 8 values of 00 
are sufficient to capture the range of model images. 



4.1 Sgr A* millimetre emission region 

The mm emission in Sgr A* is at the peak of the SED and the 
transition from optically thick to thin, and is primarily produced by 
thermal electrons. The emission comes from a photosphere where 
the optical depth r ~ 1, so that the temperature of the emitting 
electrons is given by the brightness temperature, Tj,: 



^ The simulations considered here have durations of Ri t. 
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Figure 4. Nearly face-on images from the 715h (left) and 90h (right) sim- 
ulations scaled as in Figure [T] The structure in images from untilted simu- 
lations such as 90h is mostly azimuthally symmetric (circularly symmetric 
when viewed face-on), while that in the tilted simulations such as 715h 
is highly non-axisymmetric due to the presence of standing shocks which 
dominate the mm emission in models of Sgr A*. In addition, asymmetries 
in untilted images travel at the orbital speed, whereas the standing shocks 
travel at the much slower Lense-Thirring precession speed. 



2vl 



: 6 X 10" 



Sly \ 



Jy J V 230 GHz 



A9 
40 fias 



K, (1) 



where F^g Qvo) the observed flux (specific intensity) at fre- 
quency vq, and A9 is the observed linear angular size of the emit- 
ting region sc aled to the val ue found recently from mm-VLBI ob- 
servations dPoeleman et al.ir2008l: jpish et alJl20lTh . This is effec- 
tively a lower limit for the models, since in the tilted disc mod- 
els the accretion flow can be completely optically thin at 1.3mm. 
However, the electron temperatures in these cases are still within 
a factor of 2 of Tt- This is why even for a wide range of simula- 
tions, Ti/Ts, and M , the electron tem perature is well constrained 
for both the untilted JPexter et al and tilted (ijs} models. 

In untilted simulations, the accretion flow effectively consists 
of one population of electrons - those in the inner radii near the 
midplane, where temperatures, field strengths, and particle densi- 
ties are highest. The best-fitting value of Ti /Te in untilted simula- 
tions is then that for which this electron population has Te ~ Tt- 
At higher values of Ti/Te, a larger accretion rate is required to 
match the observed flux (jV oc M'^T^ / M'^). The absorption co- 
efficient is less strongly dependent on Te than the emission coef- 
ficient (q^ oc jv/B,j oc jv/Ts in LTE), so that absorption be- 
comes increasingly important for higher Ti/T^. The accretion flow 
then develops a photosphere which grows outwards in radius with 
increasing Ti/T^. This effect is shown in the top row of Figure 
[T]for the 90h simulation. Other untilted simulations behave simi- 
la rly, except that in c onserv ative simulations (e.g., MBD o r MBQ 
in lDexteretal]|2010l from iMcKinnev & Blandfor3l2009l) the Ti 
values are larger and hence so are the favored values of Ti/T^. 

The brightness temperature argument holds for the tilted simu- 
lations as well (see Sj5j. However, heating from the standing shocks 
effectively leads to additional electron populations at varying den- 
sities, field strengths, and ion temperatures, depending on latitude. 
The densest material with the lowest ion temperature is near the 
(tilted) midplane, and is brightest for small Ti/T^. The densities 
and field strengths decrease away from the midplane, and these 
populations of electrons become more important at higher Ti/T^. 
The material closer to the midplane becomes optically thick as 
Ti/Tf, increases, but not nearly as much as in untilted simulations. 



This is because the emissivity of hotter electrons also increases, and 
these electrons are largely unobscured by optically thick material. 

Changing Ti/Te then can either lead to a growing photo- 
sphere, or it can simply increase the emissivity of electrons that 
previously had Te > Tt, causing structural changes without the 
image necessarily becoming optically thick. This effect can lead 
to a complicated dependence of the image on Ti/T^, as shown in 
the bottom two rows of Figure[T|for the 3 1 5h and 5 1 5h simulations. 
While the 3 15h image becomes increasingly optically thick with in- 
creasing Ti/Te, it does so much more slowly than the untilted sim- 
ulations. The 515h simulation begins to grow a photosphere until 
hotter electrons dominate the image at T/Te > 3, after which the 
photosphere begins to recede. These images remain qualitatively 
similar from T/Te ^ 20 - 200. 



5 SGR A* PARAMETER CONSTRAINTS 

Even with the limited coverage and sensitivity of current mm- 
VLBI data, it is possible to estimate parameters of both the black 
hole and accretion flow in S gr A* in the context of either RIAF 
teroderick et al.l l2009l 1201 ih or GRMHD models jPexter et alj 
|201(]| : IPexteJ 2()lll) . The inclination and position angles are par- 
ticularly well constrained, and are in excellent agreement between 
the two types of models. Both models assume alignment between 
the disc and black hole, which is almost certainly not the case for 
Sgr A*. 

For the tilted disc simulations with a misalignment of 15°, 
there are many possible models which provide excellent fits to cur- 
rent mm observations (reduced from a single epoch of mm- 
VLBI of < 0.6). The wide range of images is mostly due to the 
strong changes with observer time and azimuth, both of which 
have minimal effects on the structure of untilted disc images. As 
a result, the parameters of the model are essentially unconstrained 
with current observations (Figures |2] and |3]- We can formally es- 



4+1 X 10"K, all to 



timate the parameter values as: i — 80°l42o, ^ 
M = leSiiet X lO"'^M0yr-\ and T 
90 per cent confidence. All values of (/)o are allowed at this con- 
fidence. More data and a more complete set of simulations will 
be required to meaningfully constrain the parameters. All previ- 
ous models have assumed a disc tilt of zero, and have found much 
stronger constraints on Sgr A* parameters than we find for a tilt of 
15°. Previously reported parameter estimates are then invalid ex- 
cept for either very low values of tilt or spin (a < 0.3 or /3 < 15° 
are the upper limits from these simulations). 

As in untilted models, face-on inclinations are disfavored, pri- 
marily because th ey are too optically thin to explain the observed 
spectral indices (Mosc ibrodzka et al]|2009l) . A wide range of posi- 
tion angles are viable, but the location of the broad peak in prob- 
abil ity density corresponds to the range foun d for untilted mod- 
els jPexter et al.ll201(il : iBroderick et al.ll201 ll) . In each model, an 
observer azimuth is favored where velocities downstream of one 
of the non-axisymmetric shocks are approaching the observer. For 
915h, both such orientations are favored, while in the other simu- 
lations the probability distribution is asymmetric with one standing 
shock favored. 

As discussed above it is non-trivial to find all viable val- 
ues of T/Te for the tilted simulations. Since a higher accretion 
rate is required at higher T/Te to match the observed flux of Sgr 
A*, this leads to a wide range of possible accretion rates, unlike the 
untilted case. However, the electron temperature is still well con- 
strained. Lowering T/Te and increasing the accretion rate serves 
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Figure 5. Best-fitting images from each simulation (columns) with (second row) and without (top row) the effects of interstellar scattering, as well as the 
corresponding visibility amplitudes (third row) and phases (fourth row) for fo = 1.3mm (230 GHz). Image scales are the same as in Figure[T] The visibility 
amplitudes also have a dynamic range of 60, and both the visibility amplitude and phase images are 1200 X 1200AfA in size. The phase colours range 
from blue (negative) to green (zero) to red (positive). The triangles show the locations of existing mm-VLBI observations. Once convolved with interstellar 
scattering, the best-fitting images are crescents with significant but non-dominant Doppler beaming. The black hole shadow appears as a local minimum in the 
visibility amplitude, and a region of rapidly varying phase. 



Table 2. Best Fit Model Parameters 



Name 


Spin 


]\f(10-"iVlQyr-^) 


i 


« 


00 




315h 


0.3 


130 


70° 


-13° 


ir/4 


3 


5I5h 


0.5 


18 


40° 


-81° 


-n/4 


1 


715h 


0.7 


10 


90° 


-70° 


-n/2 


1 


915h 


0.9 


6 


80° 


-38° 


37r/4 


1 



to pick out different electrons, but does not change their tempera- 
ture by more than a factor of two. 



6 MODEL PROPERTIES 

Despite the lack of robust or meaningful parameter constraints, 
there are many generic properties of the images, spectra, and hght 
curves of tilted accretion disc models with important observational 
implications for Sgr A*. The best-fitting models from each simula- 
tion are representative of the range of viable possibilities from our 



fitting, and in this section these models are studied in detail. Their 
parameters are given in Table |2] 



6.1 Images 

The importance of non-axisymmetric standing shocks for Sgr A* 
images is shown in Figure|4] Images from a tilted simulation (715h, 
left panel) are dominated by material heated by the standing shocks, 
which appear as non-axisymmetric, roughly m = 2 structures. Any 
non-axisymmetric structure in unfilled images (right panel of Fig- 
ure |4j travels at the orbital speed, while the standing shocks move 
at the much slower precession speed. In addition to the choice of 
Ti/Te described above, images from tilted simulations are then pri- 
marily set by observer azimuth, which determines the viewing ori- 
entation of the standing shocks. 

The asymmetry of the tilted disc images leads to a wide range 
of complex image morphologies, which depending on the model 
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Figure 6. Visibility amplitude vs. baseline length for the best-fitting models 
interpolated to the uv locations of current and future baselines. The local 
minimum in the visibility amplitudes for 515h and 715h near 3000MA is a 
signature of the black hole shadow. 



and observer time can be dominated either by the post-shock fluid, 
dense material near the midplane of the tiUed disc, or other hot elec- 
trons. This is the primary reason that the model parameters cannot 
be well constrained thus far. Despite the additional complexity, vi- 
able images tend to be similar. This is demonstrated using the best- 
fitting models from all simulations, whose 1.3mm images (top row 
of Figure|5} are all "c rescents," as found previously for images from 
untilted simulations ( lDexteretal1l2010h . The crescent morphology 
results from a combination of strong Doppler beaming from orbital 
disc motion, which causes asymmetry between approaching and 
receding material, and gravitational lensing, which causes the back 
of the accretion flow to appear above and below the black hole in 
the image. The crescent morphology is especially apparent when 
the images are convolved with interstellar scattering (second row 
of Figure |5), which at 1.3mm blurs most of the small scale image 
structure. The black hole shadow is apparent as a central minimum 
in all images, at the orientation corresponding to light coming from 
the back of the accretion flow, which is bent above and below the 
black hole. The shadow appears in the visibility amplitude (third 
row of Figure |5} as a local minimum on the same orientation at a 
baseline whose length corresponds to the size scale of the shadow 
in the image. It also appears at the same location as a region of 
rapidly varying visibility phase (fourth row of Figure |5]. 

We can make predictions for future mm-VLBI observations by 
interpolating visibilities to the locations of future baselines between 
ALMA/ APEX in Chile, LMT in Mexico, and current telescopes: 
SMTO in Arizona, JCMT/SMA in Hawaii, and CARMA in Cali- 
fornia. The results are shown in Figure |6] for the visibility ampli- 
tudes for the best-fitting models. The models make a wide variety 
of predictions for future baselines, and the black hole shadow ap- 
pears as a local minimum in the visibility amplitude near 3000Af A 
in the 515h and 715h models. 

In addition to visibility amplitudes, phases also cany impor- 



tant information about the image structure. The closure phase is 
formed by summing phases over a triangle of baselines, and is more 
robust to instrumental errors than individual phase measurements. 
For favored (and most) position angles for all best-fitting models, 
the closure p hase is near zero, in agreement with recent mm-VLBI 
observations jFish et al.ll20T ij). In addition, the black hole shadow 
shows up as a highly variable closure phase either vs. position an- 
gle or time for position angles where the shadow is accessible. 
For nearly co-linear baselines such as those including current mm- 
VLBI telescopes, the closure phase tends to be near zero except for 
viewing orientations corresponding to the black hole shadow. How- 
ever, phase information is more interesting for other baseline trian- 
gles. Depending on viewing orientation and model, a wide range 
of possible closure phases are possible, many of which differ sig- 
nificantly from zero. Measuring closure phase in future mm-VLBI 
campaigns will therefore provide significant additional constraints 
on the models. 

At 0.87mm (Figure|7}, the other potential wavelength for fu- 
ture mm-VLBI observations, the effects of interstellar scattering 
are greatly reduced. These images have intrinsically steeper emis- 
sivities, however, so that the emission is more concentrated and 
the black hole shadow is less apparent. For these reasons, shorter 
wavelengths may be better suited for constraining parameters and 
testing accretion models than for detecting the black hole shadow. 
The structure and total flux of the accretion flow at this wavelength 
are poorly constrained, and the different simulations predict a wide 
range of both possible image structures and fluxes. For this reason, 
0.87mm VLBI may be better for constraining the models. Even 
a simultaneous spectral constraint at 0.87mm along with 1.3mm 
could distinguish between many of the possibilities. 



6. 1. 1 Implications for mm-VLBI observations 

While a crescent morphology is often found for mm tilted disc im- 
ages of Sgr A* (especially at 1.3mm), similar to previous mod- 
els of untilted discs, their time-dependent structures are completely 
different. Untilted disc images are basically static, as shown in 
the bottom row of Figure [8] Four images from the 90h simulation 
are shown spanning the simulation duration (~ 25 hours for Sgr 
A*). Even as the flux varies by ~ 50 per cent, the image struc- 
ture remains virtually unchanged. Tilted disc images, on the other 
hand, can undergo major structural changes even on relatively short 
timescales (hours, top row of Figure |8). Thus far, no significant 
structural changes have been detected with mm-VLBI, consistent 
with the untilted models. With better coverage, higher sensitivity, 
and more epochs of data, time variations in structure will become a 
powerful means to distinguish between untilted and tilted accretion 
disc models. 

Lense-Thirring precession (!j2j is another source of possible 
inter-epoch variability for tilted disc models of Sgr A*. The im- 
age structure depends strongly on the observer azimuth (Figure|9]l, 
which would effectively change as the disc precesses. Future mm- 
VLBI observations could then find evidence for this general rel- 
ativistic effect in the Galactic centre. Unfortunately, the relevant 
timescale is highly uncertain. It depends sensitively on the effec- 
tive "outer" radius of the accretion flow, its radial density profile, 
and the black hole spin. The outer radius is especially problematic, 
since both its location and interpretation are uncertain for the col- 
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Figure 7. Scatter-broadened best-fitting images at 1.3mm (top) and 0.87mm (bottom). The scaling is tlie same as in Figure [T] At 0.87mm, the images are 
easier to distinguish and the black hole shadow is less prominent. 



lisionless plasma around Sgr A*. In terms of these parameters, the 
precession time is given by tFragile et al..2007.) : 



_ 1.3 / JWt \2 



-V 

a VlOOOM/ 



(2) 



where rout is the effective outer radius. We have assumed p oc 
r~^^^ as expected for Bondi accretion and as found by hy- 
drodynamical simu lations of large-scale accretion onto Sgr A* 
dCuadra et al]|2006h . The precession time is shown in Figure [TO] as 
a function of outer radius for both the Bondi and RIAF solutions. In 
both cases, for large ranges of outer radius, the precession timescale 
is likely to be longer than the length of individual mm-VLBI cam- 
paigns (days, shaded gray region in Figure fTOt but shorter than the 
length of the entire Event Horizon Telescope experiment (years). 
This is the case, for example, for both density profiles if the ef- 
fective outer radius corresponds to the c ircular ization radius found 
for the accretion flow by Cuad ra et al.i (l2006h (orange shaded re- 
gion in Figure [TOt. Alterna tively, if the outer radius is extremely 
small (r,nax < lOOM), as IPang et alj ( 1201 ih found using MHD 
simulations with a strong large-scale magnetic field, then it could 
be possible for individual mm-VLBI epochs to span a precession 
time. 

Current data are insufficient to meaningfully test for changes 
in observer azimuth (and precession) between or within epochs. 
The prospects for detecting Lense-ThiiTing precession with Event 
Horizon Telescope observations will be considered in more detail 
in future work. 



6.2 Spectra 

Spectra from the best-fitting tilted disc models are shown in Figure 
nn The black lines show the average over observer time, while 
the dark gray envelope indicates the range. These models can- 
not explain the radio emission outside the sub-mm bump due to 
their limited spatial extent and lack of non-thermal emission. The 
IR tilted disc spectra, however, are much different than any from 
previous MHD simulations. With purely thermal electrons and no 
additional free parameters, all best-fitting models to the mm data 
are consistent with upper limits in the far- and mid-IR while pro- 
ducing the observed flux range in the near-IR (NIR) (although 



the 915h model marginally violates the mid-IR upper limit from 
iSchodel et al. 2011). 

The mean spectra from the 515h, 715h, and 915h models 
are in remarkable agreement with reported mean value s from 
ISchodel et alj ( 1201 ih and the median value from .Dodds-Eden et alj 
( I2OII ). Further, the NIR spectral indices (defined as oc 
j^-^NiR listed in Tabl e [3) from t hese mode ls are with in the 
range reported for Sg r A* dOhez et al.fc oOS : Gill essen et"ai1 [2006: 
iHomstein et alJuOOTT) . The spectral index from th e 515h model is 
in rem arkable agreement with the result found by iHomstein et alj 
( l2007l) . while those from the 715 h and 915h models agr ee with the 
redder reported spectral index of i Gillessen etal.l(l2006l) . The 315h 
spectrum underproduces the average NIR emission and its spectral 
index is steeper than observed. Other model spectra from these sim- 
ulations may provide excellent fits to the observed NIR emission in 
Sgr A* as well, and the NIR data could be included in the joint 
fitting in future studies. 

Spectra from untilted disc simulations have similar shapes 
to that shown for 90h in the left panel of Figure [12] (e.g., 
iMoscibrodzka et alj|200 ^. The NIR emission is on the exponen- 
tial tail of the synchrotron spectrum. For this reason, the untilted 
models underproduce the observed emission by several orders of 
magnitude and have much redder spectral indices than observed. 
Comparing the 90h spectrum with that from 515h (right panel of 
Figure [T2j shows that distinct electron populations (and/or a non- 
thermal distribution function) are crucial for simultaneously pro- 
ducing the mm to NIR spectrum. 

The NIR emission in the tilted models is produced by ex- 
tremely hot electrons (Te « 1 — 4 x 10^^ K) near the black hole 
(r ~ 3M). Heating from the standing shocks in the tilted disc 
models allow for multiple electron populations, some of which can 
be at very high temperatures. These can effectively mimic a power 
law tail to the synchrotron emission, especially since they are at 
significantly lower particle densities than the electrons producing 
the mm bump ( i]6.4t . This physical heating mimics the effect of 
adding a non-thermal tail to the electron distribution function (e.g.. 
Yuan, Ouataert & Naravan 2003). 

These hot electrons are also the reason we discard models with 
Ti/T^ > 3 for the 515h and 715h simulations, despite the fact 
that some of them can fit the mm-VLBI and spectral observations. 
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Figure 8. Images at four observer times for the 315h (top row) and 90h (bottom row) simulations including the effects of interstellar scattering. The scaling is 
the same as in Figure[T] The structure of images from untilted models such as 90h show very little time variability, while those from tilted models can change 
significantly on hour timescales. 
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Figure 9. Images at four observer azimuths for the 3I5h simulation including the effects of interstellar scattering. The scaling is the same as in Figure[T] The 
image structure changes significantly with observer azimuth, and these changes will be observable on the precession timescale. 




Table 3. NIR Vaiiability Parameters 



10^ 10' 
Radius (M) 

Figure 10. Precession time as a function of outer accretion flow radius 
ranging from the event horizon to the Bondi radius for radial power law 
surface density profiles corresponding to RIAF and spherical (Bondi) ac- 
cretion models. If th e effective outer acc retion flow radius is < ICp^'^M 
(gray shaded region, ICuadra et al1l2006h . then precession should occur on 
year timescales, potentially observable as changes between epochs of mm- 
VLBI observations. Only for extremely small outer radii would significant 
precession occur within a single epoch of observations (day timescales, yel- 
low shaded region). 





"NIR 


cxp /I (mJy) 




a 


315h 


-4.1 ± 1.0 


(0.7- 1.9) X 10"^ 


11 


- 18 


515h 


-0.5 ±0.1 


1.1 - 1.3 


1.8 


- 2.0 


715h 


-1.6 ±0.7 


0.5 - 0.7 


3.0 


- 3.4 


915h 


-2.2 ±0.4 


1.0 - 1.3 


1.6 


- 1.7 



Larger accretion rates are required to match the observed mm flux 
for larger Ti/T^. The peak synchrotron frequency scales as. 



while the peak flux scales roughly as 



(3) 



(4) 



so that i^c decreases for larger Ti/Te at fixed flux. This scaling 
causes the high frequency parts of the spectra to shift up in normal- 
ization and down in frequency. The approximately power law tail 
can extend into the optical/UV when Ti/Te = 1, and increasing 
its normalization leads to the violation of many IR upper limits. 
Although we do not include IR data in the fitting, we do discard 
models which violate many IR upper limits. 
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Figure 11. Best-fitting model spectra compared to Sgr A* data. The solid curves are the med ian values at each frequency over observer time, while the 
dark gray envelope shows the range. Radio data are from Falcke et al. 1 1998, open circles) and lAn et al.l 12005. filled diamonds). The models are fit to the 
mm data (filled squares, Marrone 2006), but are in good agreement with the mean/median (orange and purple, Schodel et al.ll201 11 : Dodds-Eden et al 201 ll, 
respectively) and flaring (pink, Genzel et al. 2003) NIR data with no additional free parameters. They also satisfy observational upper limits in the mid-IR 
(blue upper limits, Melia & Falcke 200 1 , a nd references there in). At lower fre quencies the non-thermal emission is produced outside of the simulation domain, 
either in the accretion flow I I Yuan. Ouataert & Naravanll2003l) or in a short jet ffalcke & Markofj2000l) . 



6.3 Millimetre to infrared variability 

Radiative models of Sgr A* based on simulations are 
thus far the onl y to self-consistently include va riabil- 
ity information ( Goldston, Ouataert & Igumenshchev '2005'; 
Dexter. Agol & Fragi le 2009; Dexter etal. 2010; Dolence et al. 
2012h . Best-fitting models to mm-VLBI and spectral infor- 
mation from untilted GRMHD simulations naturally produce 
millimeter variability qualitat i vely consistent w ith that observed 
toexter. Ag ol & Fragilel l2009t iDexter et al.ll2010) . This variability 
is due to relatively global (m = 0, 1) fluctuations in the particle 
density and magnetic field strength, and is strongly correlated 
with the accretion rate. We find similar mm light curves for the 
best-fitting tilted disc models, as shown in Figure [T3] All of the 
best-fitting light curves exhibit variability on hour timescales with 

< 50 per cent amplitudes, roughly in agreement with observations. 
Millimeter "flares" would be identified as two per light curve, 
consistent with their observed frequency in Sgr A*. However, 
tilted disc models produce larger amplitude and lon ger duration 
flares than their untilted counterparts (cf. Figure 6 of iDexter et alj 
I2OIQ) . The variability is still strongly correlated with the accretion 
rate, which also varies more in tilted disc simulations. 

As discussed above, tilted disc models naturally produce 
roughly the correct mean levels of IR emission, very different 
from previous numeri cal accretion flow models based solely on 
thermal electrons (e.g. .'Goldsto n. Ouataert & Igumenshchev|[200^ ; 
iMoscibrodzka et al. 2009). The H (l.e/xm), K (2.2^m), and L 
(3.6/im) band light curves from the best-fitting models to the mm 
data are shown in Figure[T4] In all cases, the NIR emission is highly 
variable, with factor of ~ 5 — 20 increases over the mean level on 

< Ihr timescales, consistent with observations. The time sampling 
in the light curves is 1 — 10 minutes, depending on the simulation. 
The 915h simulation, with 1 minute sampling, easily resolves the 
NIR flux variations. The shortest significant flux variations (factors 
of > 2) occur on timescales ~ 5 — 15 minutes. Substructure in 
the 315h and 715h light curves, with 10 minute sampling, is only 
marginally resolved. 



The various models differ substantially in the predicted NIR 
spectral slope (see §and Table |3]l. Similarly, two of the spectral 
slopes are time variable (315h and 715h), while the others are 
nearly constant (515h and 915h). The flux distributions of all of 
the light curves are fairly well described by a log-normal: 



■ exp 



2a2 



(5) 



with parameters /i and a. The median is exp(p), and a represents 
a multiplicative standard deviation. Thi s model provides a better 
description than a standard Gaussian, as iDodds-Eden etal.l ( l201lh 
found for Sgr A* data. The parameters inferred from our K-band 
light curves are given in Table |3] The ranges represent the results 
of fitting with various bin sizes. 

Power spectra from all the tilted disc models are qualitatively 
similar to that for 915h (thick solid curve in Figure [13]. The slope is 
fairly shallow (slope ~ —1) on timescales significantly longer than 
that of the marginally stable orbit (> 10 min), altho ugh the effec- 
tive inner disc edge is well outside of this location jFragilell2009l; 
IDexter &Fragilell201ll) . On shorter timescales the power spectrum 
steepens to a slope ~ —4.5. Fit with a s ingle power law, the slope is 
^_-2.3, consistent with observations ( IMever et al ]|2008l : lDo et Zl 
I2OO9I) . 

iDolence et alj j201 j) found a broken power law power spec- 
trum in a model of Sgr A* from an untilted GRMHD simula- 
tion with a similar break frequency. However, we find no signifi- 
cant quasi-periodic oscillations similar to what they reported near 
the frequency of the marginally stable orbit. Using many cam- 
eras spaced in observer azimuth, they decomposed the power spec- 
trum into azimuthal modes (their Eq. 1 and Fig. 1). A compari- 
son is shown in Figure [TS] Unlike in their models, the m = 
mode is dominant, with m = 1, 2 only contributing on the short- 
est timescales. These discrepancies could be due to differences in 
the simulations (either numerical algorithms or duration), or due to 
physical differences between tilted and untilted accretion discs. 
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Figure 12. Spectra from best-fitting untilted (90h. left) and tilted (5 15h, right) simulations. Sgr A* data are the same as in Figure [TT] In both cases the spectra 
are fit to the green sub-mm points. Multiple electron populations from shock heating in tilted discs can naturally produce the observed NIR emission, which 
is underproduced by Ki 2 orders of magnitude in untilted simulations. Other untilted simulations exhibit larger amplitude variability than 90h, but its spectral 
shape is representative. 
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Figure 13. Light curves from the best-fitting models at 1.3mm (230 GHz, solid lines) and 0.4mm (690 GHz, dotted lines). The times of best fit to the visibility 
data are shown as upward fiiangles, while the ranges of times averaged for the best fit to the spectral data are bracketed by left and light facing arrows. The 
simulated light curves exhibit 50 per cent variability on hour timescales, similar to the untilted case and in qualitative agreement with observations. The mm 
spectral index increases with flux, because the optical depth is lower at higher frequency. 



6.4 Physical origin of time-variable mm and NIR emission 

The mm emission in the tilted disc models arises from hot, dense 
electrons in the inner radii of the simulations. Typical parameters of 
mm electrons are similar to those in untilted discs: n ~ lO^cm"^, 
B ~ 50G, Te ~ 6 X 10^°K. As discussed above, this tends to be 
where densities and field strengths are highest: either in the post- 
shock fluid or in the midplane of the inner radii of the tilted ac- 
cretion flow. Due to the asymmetry in the shocks and subsequent 
infall, the emission geometry is different from untilted simulations. 
The mm variability in the tilted disc models is, however, essen- 
tially the same as in untilted models. The variability is caused by 



relatively global fluctuations in the particle density and magnetic 
field strength. Most of the power in these fluid variables as well as 
in the emissivity is in the m = mode, although there is also sig- 
nificant power in m = 1, 2, particularly in the form of spiral waves. 
The 0.4mm (690 GHz) light curve is strongly correlated with the 
accretion rate through the inner boundary, suggesting that at least 
the m = power is from global changes in the accretion rate. The 
global nature of the variability may explain its similarities in tilted 
and untilted models, despite significant differences in their dynam- 
ics. 
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Figure 14. Near-infrared light curves from the best-fitting models in the H, K, and L bands. All of the simulated NIR light curves are highly variable on short 
(< hr timescales). The 315h simulation underproduces the observed NIR flux and its vaiiability. The 515h, 715h, and 915h models produce the observed flux 
range and exhibit "flares" with roughly the right amplitudes (factors < 30) and timescales (~ Ihr). The 315h and 715h spectral slopes are time-variable, 
while those from the 515h and 915h models are nearly constant with NIR flux. 



iDexter et al. found that a simple approximate emissiv- 

ity, ji, oc nB^, captured the shape of the mm light curves from un- 
filled simulations of Sgr A*. This is the expected form for isother- 
mal synchrotron emission with Te ~ 6 x 10^" K, and demonstrates 
that temperature fluctuations are not important for driving the vari- 
ability. We find the same result for the tilted disc models. This vari- 
ability mechanism is not magnetic reconnection, since the tilted 
simulations do not capture heat lost by grid-scale magnetic dissipa- 
tion. It is also not orbiting hot spots, as the m = power dominates 
the fluctuations. 

The NIR emission in the models is entirely optically thin, 
and dominated by very small regions at small radius (r ~ 3M) 
in the simulation with lower densities (n ~ lO^cm"'^), simi- 
lar magnetic field strengths {B ~ 50G), and hot temperatures 
(Te > This location is fairly insensitive to spin, and is 

inside the marginally stable orbit for all but the 915h simulation. 
These are lower particle densities and higher electron temperatures 
than previously reported for IR-emitti ng electrons from axisym - 
metric, untilted GRMHD simulations jMoscibrodzka et al.|[2009l) . 
In their simulations, the NIR emission was dominated by heating 
from numerical magnetic reconnection in current sheets, which is 
thought to be an artifact of axisymmetry. In our model, the NIR- 
emitting electrons are heated by the non-axisymmetric standing 
shocks. This is the only source of entropy generation in these sim- 
ulations, since heat from numerical reconnection is lost from the 
grid. For example, the maximum temperature obtained anywhere 
in the 90h simulation is < lO^^K. 

The NIR "flares" in our models are due to fluctuations in the 
particle density, field strength, and electron temperature of these 
hot post-shock electrons. The 915h K-band light curve can be re- 
covered approximately with an emissivity j„ oc nB^, as long as 
only hot electrons are included (Te > 5 x 10^^ K). In the 515h 
and 715h simulations (with hotter electrons), the emissivity can be 
approximated as oc nBT^. Temperature fluctuations are impor- 
tant in the 515h and 715h simulations, but not in 915h. The ap- 
proximate forms can be understood by the variations in emissivity 
expected from fluctuations in the individual quantities at the appro- 
priate ratio of u/vc- This ratio is near 1 for the hotter 515h and 
715h simulations, and ~ 10 — 20 for 915h. Closer \.o = Vc, the 
emissivity is less strongly dependent on B and Te . 

General relativistic effects are also important for the NIR light 
curves from 515h and 715h. Many of the flares in these models 
correspond to extreme gravitational lensing events, where emission 



from hot electrons behind the black hole is lensed into a bright ring. 
An example from the 5 15h model is shown in Figure[T6] A flare oc- 
curs between the first and second panels when hot gas behind the 
black hole is strongly lensed, concentrating its emission in a bright 
ring at the circular photon orbit. In the subsequent frames, this hot 
material passes by the black hole and moves to larger radius, lead- 
ing to a decay of the flare. The lensing acts as a filter, enhancing 
the emission from electrons at larger radius (r > 3Af ), and can 
produce modulations in the light curves of order the observed flux. 

The NIR image structure is highly variable, and so is its cen- 
troid. The K-band centroid position is shown in Figure[T7] The typ- 
ical excursions are ~ 30 — 50/ias in all best-fitting models except 
that from the 915h simulation. The position wander is only weakly 
correlated with the NIR flux. All periods of rapid variability and 
very high NIR fluxes also exhibit significant centroid motions, but 
not all periods of large position wander correspond to NIR flares. 

The NIR spectral index is set by the critical frequency for syn- 
chrotron emission. 



This primarily depends on electron temperature, which varies be- 
tween simulations. The NIR electrons are coldest in the 315h sim- 
ulation, and hottest in the 515h simulation. Somewhat surprisingly, 
the shock heating is less effective in the 715h and 915h simulations. 
lDexter&Fragii3(i2011) argued that the non-axisymmetric stand- 
ing shocks are stronger in higher spin simulations. Colder NIR elec- 
trons and lower critical frequencies lead to more negative spectral 
indices. These also tend to be more variable, since the emissivity 
depends strongly on frequency when v jvc 3> 1. When v jvc — 1, 
as for the 515h simulation, the emissivity is fairly flat, and the spec- 
tral index is nearly constant. 

The mm and NIR light curves are uncorrelated in this model, 
since they arise from different populations of electrons. The full 
690 GHz and K-band light curves and their cross-correlation are 
shown in Figure [18] Although the full duration light curves are 
uncorrelated, 6 hour segments can produce spurious peaks in the 
cross-co rrelation w ith hour lags, similar to those reported for Sgr 
A* (e.g., lYusef-Zadeh et aL2008,) . 
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Figure 15. NIR K-band (2.2/im) power spectrum (thick solid curve) from 
the best-fitting 915h model, averaged over 10 evenly-spaced light curve 
segments. The power spectrum is further decomposed in azimuthal mode, 
as done in Dolence et al. ( 2012). Both the overall power spectrum and those 
from individual azimtuhal modes are well described by broken power laws 



with break frequencies 
~ —1 and —4.5. 



7 DISCUSSION 



lOmin . The full power spectrum has slopes 



Geometrically thick, low luminosity accretion flows do not align 
with their central black holes. The current source of accretion onto 
Sgr A* is winds from massive stars, whose net angular momentum 
is almost certainly unrelated to the spin axis of the black hole. Both 
the current accretion rates onto Sgr A* at the Bondi radius and in 
the inner radii are orders of magnitude smaller than required for 
sufficient mass accretion to have changed the black hole spin ori- 
entation over a Hubble time. Therefore, it is highly unlikely that 
the black hole spin axis in Sgr A* aligns with the angular momen- 
tum of the accreting gas. We have constructed the first radiative 
models of Sgr A* from misaligned (tilted) black hole accretion 
flows by performing radiative transfer calculations on data from 
GRMHD simulations (Table [T]). These models provide an excel- 
lent description of existing observations, and are significantly dif- 
ferent from any theoretical model previously considered. Even for 
the modest 15° misaligrmient considered between black hole spin 
and accretion angular momentum axes, the (thermo)dynamics of 
the accretion flow change drastically. Non-axisymmetric shocks 
from the convergence of eccentric orbits in the warped disc give 
rise to asymmetric image morphologies, leading to the breakdown 
of many simplifying assumptions in the untilted models. The ob- 
server azimuthal viewing angle becomes an important new param- 
eter, and the dependence of observables on both the viewing ge- 
ometry and the ion-electron temperature ratio become non-trivial. 
The combination of these factors greatly expands the possible pa- 
rameter space for tilted disc models relative to those assumed in all 
previous models of Sgr A*. 

There are numerous implications. The parameters of both 



the black hole and the viewing geometry are poorly constrained 
(Figures |2] and [3} by mm-VLBI and spectral observations. Pre- 
viously reported parameter cons traints (^Br oderick et al. 2009, 
201ll; iM oscibrodzka et al. 2009 ; IPexter, Agol & Fragile 2009; 



Dexter et al.i i2010i ; iShcherbakov. Penna & McKinnevi i201Q ) only 



apply if the black hole spin or tilt values are extremely small. Upper 
limits from this work are a < 0.3 or /? < 15°, although qualita- 
tively similar behavi or has been seen in a simulation with jS = 10° 
( iHenisev et ai]|2009l) . The best-fitting models (Table[2ll to the mm 
data sort into two groups at low (10~^ — lO~*M0yr~^) and high 

(> 10~^MQyr^^) time-averaged mass accretion rates^ 

The "crescent" black hole images reported by ' Dexter et al.l 

1 I2OIOI) are also found from these best-fitting tilted images at 
1.3mm, especially when including the effects of interstellar scat- 
tering. For this reason, the black hole shadow is still expected 
to be accessible, although huge uncertainties in the viewing ge- 
ometry for these models prevent a robust prediction for baselines 
where it should be seen. In two of the best-fitting models (see 
Figure |6l, the shadow is visible on baselines between telescopes 
in Chile (ALMA/APEX) and Mexico (LMT), but a wide vari- 
ety of other geometries are nearly as likely. The mm variabil- 
ity properties repo rted previously dDexter, Agol & Fragilel |200^; 
iDexter et al.l l201(ii) qualitatively hold for the tilted disc models, 
although their variability amplitudes and characteristic timescales 
seem somewhat larger (Figure [13) • 

Although images from tilted discs are crescents at 1.3mm, 
at 0.87mm the images can become double-peaked due to emis- 
sion from near the standing shocks (915h model in Figure |7). 
The structure of the tilted images is also highly time-variable, 
in stark contrast to the static morphology in the untilted case 
(Figure [§}. Tilted and untilted models can then be distinguished 
either with higher frequency mm-VLBl or simply with many 
epochs of observations using the current array. Polarized mm- 
VLBI observations may also help to dis tinguish between models. 
IShcherbakov, Perma & McKirmevI I2OIQ ) found that spectropolari- 
metric data contain comparable constraining power as current mm- 
VLBI observations. 

Perhaps the most important difference is that the low accre- 
tion rate best-fitting models (those from the 515h, 715h, and 915h 
simulations) naturally reproduce IR spectral (Figure II 11 and vari- 
ability (Figure [14) properties consistent with observations, while 
those at high mass accretion rates (from the 315h simulation) un- 
derproduce the observed flux. All m odels satisfy the observed up- 
per limits on mid-infrared emission ( Telesco. Davidson & Wemen 
1 19961; ICoteraetalJll999l ; ISchodel et al.ll201lh ! All three low accre- 
tion rate models are in excellent agreement with the reported mean 
NIR fluxes from Sgr A* (Schod el etaLll201ll ; iDodds-Eden et alj 
[201 Ij). Flux distributions from the 515h and 915h models are well 
described by a log-normal, and the inferred pa rameters (Table [3) 
are in excellent agreement with those found by iDodds-Eden et alj 
( I20T J). Both the mean and time-dependence of the NIR spectral in- 
dex varies considerably between models. The 515h spectral index 
is constant in time with a va lue of onir = — .5± 0.1, nearly iden- 
tical to the observations of iGhez et al.l ( l2005h and iHomstein et al] 
1 I2OO7I) . The 715h and 915h light curves have much redder spec- 
tral indices, and the 715h spectral index increases with flux. These 
values are in good agreement with Gillessen et al. (2006). Because 
of the short duration of our light curves (~ 25hr), it is unclear 
whether the 515h, 715h, and 915h light curves are best interpreted 
as different states of the same theoretical model, or as a means to 
distinguish between the models. 

In any event, this is the first model of Sgr A* that self- 
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Figure 16. H-band (1.6/xm) images from the 515h model spaced ~ 5min apart during an NIR "flai'e" with a factor of ~ 4 flux variation between the first and 
second images. The images are scaled logarithmically with a dynamic range of 2048, and the panel size is 85 X 85/ias. The flare emission is clearly lensed 
and concentrated at the photon ring, coiTesponding to the unstable photon orbit in the Kerr spacetime. 
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Figure 17. X (solid) and Y (dotted) centroid positions for the K-band 
(2.2/im) images. Periods of rapid centroid movement of ~ 30 — 50/tas 
occur every few hours in all models but 915h, and should be detectable by 
the future VLT instrument GRAVITY. 



consistently produces the observed time- variable mm to NIR emis- 
sion. The two are completely uncorrelated. The mm variations 
are due to magnetic turbulence and are strongly correlated with 
the accretion rate. The NIR variations are caused by particle den- 
sity, field strength, and temperature variations in tenuous, hot gas 
heated by the non-axisymmetric standing shocks, and can be signif- 
icantly modulated by extreme gravitational lensing. Short segments 
of these light curves can, however, still lead to the reported weak 
correlations (Figure[T8). We caution therefore that cross- correlation 
coefficients of < 0.6 between wavebands in single epoch observa- 
tions do not necessarily imply a common physical origin. 

This NIR emission model is completely novel. It is not 
due to orbiting or e xpanding "hotspots" ([B roderick & Loeb 2005 



I 



lEckart et al.l bOOd lYusef-Zadeh et al.l I2OO& : .Eckart et al.. ,201 
However, significant structural changes can still be seen during 
large changes in NIR flux (Figure [TTll. It is similar in spirit to mod- 
els that invoke magnetic reconnection events tPodds-Eden et alj 
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Figure 18. Top: Light curves from the 515h simulation at 0.4mm (red) 
and 3.6^m (blue). Bottom: Cross-correlation of the full light curves (solid 
curve) and four evenly spaced segments (dashed curves). For these short 
(several hour) durations, lags can be found at correlations of ~ 0.4 — 0.6 
as reported in the literature (e.g., lYusef-Zadeh et al.lbooil) . despite the fact 
that the full light curves are uncorrelated. 



I2OIOI) to accelerate electrons to the observed high energies. But 
magnetic reconnection cannot heat particles in these simulations. 
Instead, the heating is from the non-axisymmetric standing shocks 
in the tilted accretion flow. This emission is also completely dif- 
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ferent from that described by iDolence et alj ( 1201 2b . who detected 
quasi-periodic oscillations in NIR/X-ray light curves of Sgr A* 
from an untilted GRMHD simulation. Their model underproduces 
the observed NIR flux by an order of magnitude, simil ar to previous 
models from untilted simulations jPexter et al.ll20ld) . In addition, 
the NIR/X-ray variability in their model is strongly correlated with 
that in the mm, in conflict with observations. 

These radiative calculations assume that cooling is negligible 
for the fluid (thermo)dynamics. However, the cooling timescale for 
the NIR-emitting electrons is short: 

At r ~ 3A/, where the emission is produced, this timescale is still 
larger than the dynamical time (tdyn — 2 min). If this is where 
the particle heating occurs, ignoring the effects of cooling on the 
electron temperature may be a safe approximation. We have also 
assumed a thermal electron distribution function, which is likely in- 
correct. Both of these limitations can be addressed in future calcula- 
tions by following particle trajectories through the standing shocks 
and self-consistently accounting for heating and cooling. 

We have only included synchrotron radiation. Including in- 
verse Compton sca ttering via Monte Carlo calculations (e.g., 
iDoIence et al.l [20l3) would allow the calculation of X-ray light 
curves, and in particular the X-ray luminosity during NIR flares. 
We can crudely estimate this using typical values of fluid variables 
during rapid changes in NIR flux: Te = 2 x 10^^ K, R = IM 
and n = lO'^'cm"^ (low, H-band flux of 1 mjy) to = 3M and 
n = 6 X lO'^cm"^ (high, H-band flux of 10 mJy). For this tem- 
perature, the Compton peak is at ~ IkeV. The ratios of Compton 
to sy nchrotron bolomet ric luminosities estimated using the formu- 
lae in lEsin et alj ( Il996h are Lc /Ls ~ 0.1 (low) and Lc /Ls ~ 2 
(high). At those H-band fluxes, the estimated X-ray luminosities 
are: Lx ~ lO^^erg s"^ (low) and Lx ~ 10^^ erg (high). 
These num bers are in rough agreement with observations of X-ray 
flares (e.g.. lBaganoff et alll200lLl2003b . 

The effects of radiative cooling on multi-wavelength spectra 
from axisymmetric (aligned) GRMHD sim ulations have be e n stud - 
ied by Dibi et al. (2012, submitted) and IPrappeau et al] ( l20Ill) . 
They find that adding radiation in post-processing, as done here, 
is valid for accretion rates M < lO~^M0yr~^. All of our best- 
fitting models satisfy this criterion (Table |2) except for that from 
315h, where the accretion rate is slightly larger than this value. 
Radiative cooling may be important for this and similar models, 
although it is unknown whether this accretion rate limit applies to 
tilted discs. 

We have only considered four simulations, all with a tilt an- 
gle of (3 = 15°. The non-conservative GRMHD algorithm used in 
these simulations leads to a scale height H/r ~ 0.1, whereas the 
accretion flow onto Sgr A* may have H/r ~ 1. We are unable 
to probe the dependence of our results on either tilt angle or scale 
height, both of which may change the formation and strength of the 
standing shocks. It is unclear whether the models considered here 
are representative of the range of possibilities from tilted black hole 
accretion discs. 

All of these simulations failed to conserve total energy. We 
previously found that magnetic energy lost to grid-scale numeri- 
cal reconnection leads to a radiatively efficien t effective luminos- 
ity, ~ O.lMc^ jPexter. Agol & Fragilell2009l) . Although untilted 
models based on these simulations give very similar re sults to those 
that recapture this energy as heat dPexter et al.ll2010h . it is unclear 



whether that is the case for tilted discs. The lack of energy conser- 
vation means that the only form of entropy generation is through 
shock heating. Numerical dissipation of magnetic fields would pro- 
vide another heat source, and could potentially lessen the dramatic 
effects seen here that are caused by the standing shocks. This issue 
can be addressed by future energy-conserving tilted simulations. 

If standing shocks forming from eccentric orbits in a tilted, 
warped disc provide the heating mechanism for the observed NIR 
emission, then the inner disc must be tilted with respect to the black 
hole spin axis. Then the Galactic centre black hole must have non- 
zero spin. Furthermore, recent general relativistic, semi-analytic 
models of geometrically thick tilted accretion discs find that the in- 
ner d isc aligns with the black hole for a < jZhuravlev & Ivanovl 
l20Ilh . Therefore, if this proposed model for the NIR emission is 
correct, the black hole in Sgr A* must have positive spin (a > 0). 



8 SUMMARY 

There is no reason to expect alignment between the black hole spin 
axis and the angular momentum axis of accreting gas in radiatively 
inefficient sources. We have calculated time-variable mm to NIR 
emission from the only published numerical models of misaligned 
("tilted") black hole accretion discs and compared them to obser- 
vations of Sgr A*. Our main results are: 

(i) Tilted disc models provide an excellent description of exist- 
ing mm-VLBI and spectral observations of Sgr A*. 

(ii) Previously reported parameter constraints from fitting mod- 
els to mm-VLBI, spectral, and/or polarization measurements are 
only valid if the spin or tilt value is very small (upper limits of 
a < 0.3,13 < 15°). 

(iii) Predicted images for mm-VLBI are still crescents as in the 
aligned case, and the black hole shadow, direct evidence for the 
existence of an event horizon, is still potentially visible on baselines 
between Mexico and Chile. 

(iv) Heating from non-axisymmetric standing shocks can natu- 
rally produce the observed NIR emission. This emission is uncoiTe- 
lated with the mm variability, which also agrees with observations 
and is strongly correlated with the mass accretion rate. The X-ray 
flares in this picture arise from Compton upscattering of NIR seed 
photons. 

(v) We predict several observational signatures of disc tilt: a) 
structural changes observable with mm-VLBI from strong asym- 
metry in the emission region on hour (dynamical) and/or day- year 
(Lense-Thirring precession) timescales; b) bimodal structures in 
345 GHz mm-VLBI images from the presence of standing shocks 
in the accretion flow; and c) NIR centroid motions of 30 — 50/ias 
from extreme gravitational lensing during flares. 

(vi) If the NIR emission is caused by heating from standing 
shocks in a tilted accretion flow, then the spin parameter of Sgr 
A* must satisfy a > 0. 
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